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Abstract 


The  method  of  determining  the  thickness  of  diffuse  phase 
boundary  with  density  profile  governed  by  equilibrium  condition?, 
is  proposed.  It  follows  the  well-known  procedure  of  analyzing  the 
deviations  from  Porod's  law.  Errors  in  the  obtained  boundary  thick¬ 
ness,  due  to  the  statistical  scatter  in  the  intensity  data  and  to 
the  difficulty  of  separating  the  effect  of  density  fluctuations  with¬ 
in  the  phases,  are  examined.  For  this  purpose,  scattering  curves 
are  synthesized  on  the  basis  of  a  well-defined  model  structure  with 
known  boundary  thicknesses.  These  synthesized  curves,  when  analyzed 
according  to  the  proposed  method,  yield  the  correct  boundary  thick¬ 
ness  under  favorable  conditions,  but  are  also  shown  to  lead  to  very 
erroneous  results  in  some  cases. 


I.  Introduction 


The  intensity  of  small-angle  X-ray  scattering  (SAXS)  fcv  iso- 
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tropic  materials  having  ideal  two  phase  structure  decreases  as  s 
at  large  s  Cs  =  2sin0/l] ,  according  to  the  theory  by  Porod  (1?51, 
1952a,  b) .  Deviations  from  this  Porod  law  can  occur  in  practice 
since  the  boundaries  between  the  phases  may  not  be  perfectly  sharp, 
and  the  electron  densities  within  the  phases  contain  local  fluctu¬ 
ations.  The  potential  utility  of  determining  the  phase  boundary 
thickness,  by  examination  of  the  deviation  of  observed  SAXS  inten¬ 
sity  from  Porod 's  law,  was  first  pointed  out  by  P.uland  (1971)  and 
was  since  examined  by  a  number  of  workers  (Vonk  1973;  Ruland  1974; 
Hashimoto,  Todo ,  Itoi  and  Kawai  1977;  Todo,  Hashixnoto  and  Kawai 
1978;  Hashimoto,  Shibayama  and  Kawai  1980;  Koberstein,  Morra  and 
Stein  1980;  Roe,  Fishkis  and  Chang  1931). 

According  to  Ruland  (1971),  the  observed  intensity  I0bs  ^ 
can  be  expressed  as 

IobsCs)  =  Iid(s)H2Cs)  +  Ib  (1) 

where  1^(5)  is  the  intensity  due  to  ideal  two  phase  structure, 

H  (s)  is  the  factor  arising  from  diffuseness  of  the  boundaries, 
and  Ij.  is  the  'background'  intensity  due  to  local  electron-density 
fluctuation  within  the  phases.  If  the  electron-density  profile 
across  the  boundary  approximates  a  convolution  product  of  a  step 
function  and  a  gaussian  function,  H  (s)  takes  the  form 

K2  (s)  =  exp  (-4t2c^s^ )  (gaussian)  (2) 

where  the  variance  a  of  the  gaussian  function  embodies  the  diffuse¬ 
ness  of  the  boundary.  Statistical-mechanical  studies  of  inhomo¬ 
geneous  liquid  systems  (Cafcn  and  Hilliard  1953)  ,  polymer  mixtures 
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r 


CHelfand  and  Tagami  1971;  Hopper  and  Uhlmann  1974;  Roe  1975)  and 
block  copolymers  CHelfand  1975)  showed  that  the  density  profile 
across  the  phase  boundary  attained  under  equilibrium  conditions 
is  given  by 

gCxj  -  tanhCiTx//rr  cr)  (3) 

In  this  case  the  corresponding  Cs)  function  is  (Ruland  1980;  Roe, 
Fishkis  and  Chang  1981) 

Cs)  =  CvT^iras)  ^csch^  Cj^T^TTas)  C equilibrium)  (4) 

The  boundary  thickness  tor  the  value  of  a)  is  determined  by 
analysing  the  observed  intensity  in  the  Porod's  region  according 
to  eq.  CU ,  and  for  this  purpose  a  suitable  graphical  or  numerical 
method  has  to  be  devised.  The  simplest  is  to  expand  H  ts)  in  a 
series  and  retain  the  first  two  terms  only.  Kober stein,  Morra  and 
Stein  C19801  showed  that  such  a  procedure  is  likely  to  lead  to  an 
underestimate  of  cr  when  the  analysis  is  performed  over  a  range  of 
s  in  which  truncation  of  higher  order  terms  is  not  justified.  They 
instead  proposed  a  procedure,  applicable  when  Cs)  is  given  by 
ec.  C2),  which  is  valid  even  for  s  greatly  exceeding  1/c.  One  of 
the  purposes  of  the  present  work  is  to  present  a  similar  procedure 
which  is  applicable  when  H  Cs)  is  given  by  eq.  C4)  instead  of  eq.  (2). 

The  main  practical  difficulty  in  utilizing  eq.  Cl)  for  esti¬ 
mation  of  a  arises  from  uncertainties  introduced  by  the  errors  in 
the  observed  data.  The  effect  of  diffuse  boundaries,  represented 
by  H  (s) ,  introduces  a  negative  deviation  from  Porod's  law,  while 
the  'background'  1^  contributes  a  positive  deviation.  These  two 


effects  therefore  tend  to  cancel  each,  other,  and  meaningful  sepa¬ 
ration  of  one  from  the  other  is  often  found  difficult  because  of 


the  statistical  error  in  the  data,  which,  in  the  Porod's  region  is 
usually  appreciable.  It  is  desirable  to  have  an  estimate  of  the 
error  in  the  evaluated  a  value  arising  in  practice  as  a  result  of 
the  statistical  error  in  the  data .  For  this  purpose  in  this  work 
a  series  of  'theoretical'  scattering  curves  based  on  a  model  struc¬ 
ture  are  constructed,  and  are  then  given  random  statistical  fluctu¬ 
ations  to  make  then  simulate  experimentally  determined  data.  These 
synthesized  scattering  curves  are  thereafter  analyzed  by  the  method 
proposed  for  the  determination  of  a  values,  and  the  values  thus  ob¬ 
tained  are  compared  with  the  values  initially  chosen  for  the  syn¬ 
thesized  curves. 


-3- 


II.  Procedure  Applicable  to  Equilibrium  Boundary  Profile 

When  SAXS  data  are  collected  with  instruments  having  slit 
collimation,  the  effect  of  the  slit-smearing  has  to  be  properly 
taken  into  account.  Because  of  the  possible  amplification  of  the 
error  in  the  data  in  the  process  of  desmearing,  the  procedure  for 
the  determination  of  the  boundary  thickness  should  preferably 
make  use  of  the  smeared  data  directly.  The  slit-smeared  inten¬ 
sity  Its)  is  given  by 

ICS)  =  2^  W£tu)  Iobs  (Vs"2+U2)  du  C5) 

where  W^lu)  is  the  slit-length  weighting  function.  In  what  follows, 
we  will  consider  only  the  cases  in  which  the  inr inite-slit  approxi¬ 
mation  is  valid.  tu)  is  then  a  constant  independent  of  u. 

The  intensity  of  ideal  two  phase  structure  in  the  Porod's 
region  is  given  by 

lid  ts)  -  CAp)2S/8tt3s4  (6) 

where  Ap  is  the  difference  in  the  electron  density  between  the  two 
phases  and  S  is  the  total  phase  boundary  area.  For  the  density  pro¬ 
file  across  the  phase  boundary  determined  by  equilibrium  conditions, 
H2  Cs)  is  given  by  eq.  (4),  and  therefore  the  slit-smeared  intensity 
in  the  Porod's  region  can  be  written  as 

ICs)  =  Ag  lx)  +  Ib  (7) 

where  ib  is  the  smeared  'background'  intensity,  and 

A  =  36ttSCAp)2W£  (0)a4  C8) 

x  =  /12-rras  (9) 
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2  ? 

g  U1  -  csch.  (x)/x 

(101 

g  Cx]  =  <j(/x2j-u^) 

(11) 

The  integration  indicated  in  ec.  (H)  could  not  be  evaluated 
analytically,  and  in  this  work  g  (x)  was  calculated  numerically . 

For  x  between  0.8.  and  6.0,  g(x)  is  approximated  well  by 

gCx)  3  3.3x~1,8e_1-9x  C12) 

The  root-mean-square  error  of  this  approximation  over  the  above- 

mentioned  range  of  x  is  2%  and  the  maximum  error  at  the  two  ends  of 

1  8 

the  range  is  4%.  Fig.  1  shows  the  plot  of  lnlx  *  g(x)]  against  x, 

illustrating  the  good  linearity  in  the  indicated  range  of  x. 

Evaluation  of  a  from  the  observed  intensity  I(s)  would  then 

^18 

proceed  as  follows.  A  plot  of  ln{{I(s)  -  ‘  ]  against  s  is 

constructed,  and  the  slope  of  the  linear  portion  of  the  plot  is 
determined.  The  value  of  a  is  then  given  by 

a  =  -slope/ (1. 9/17^)  C13) 

Alternatively,  one  could  evaluate  the  effective  thickness  t  by 

t  =  -slope  (2/1 . 9w2)  (14) 

where  t,  defined  as  the  reciprocal  of  the  density  gradient  at  the 
center  of  the  boundary  layer,  is  related  to  c  by 

t  =»  (4/T/t r)a,  (15) 

After  cr  has  been  evaluated,  one  may  calculate  x  by  eq.  (9) 
to  see  whether  the  approximation  (12)  is  justified  in  the  range 
of  s  at  which  the  slope  was  evaluated.  In  the  (unlikely)  event 


that  the  x  range  starts  at  a  value  smaller  than  0.3,  the  following 
alternative  approximation  may  be  employed. 


~  _  *  -2.5  -1 . 4x 

gCx}  ~  2.  Ox  e  C16) 

Eq.  Cl 6)  is  valid  for  x  between  0.3  and  2.1,  with,  the  root-mean- 
•square  error  equal  to  2%  and  the  maximum  error  equal  to  4%.  The 
required  plot  is  then  ln[{Its)  -  Ifc}s  '  ]  against  s. 

III.  Synthesized  Scattering  Curves 

It  would  be  a  difficult  task  to  derive  an  analytical  expression, 
in  general  terms,  for  the  propagation  of  error  from  the  initial  data 
to  the  final  a  value  obtained.  Instead,  we  take  a  simpler  approach 
of  synthesizing  a  model  scattering  curve  and  then  applying  the  pro¬ 
posed  procedure  to  it  to  see  how  closely  the  known  a  value  is  re¬ 
covered.  In  constructing  the  model  scattering  curves,  effort  is 
given  to  make  them  as  realistic  as  possible.  The  physical  constants 
required  for  their  calculation  are  given  values  typically  found  in 
real  systems  of  interest. 

Two  types  of  scattering  curves  are  constructed.  In  the  first, 
the  model  structure  assumed  is  one  in  which  spherical  particles  of 
one  phase  are  imbedded  in  a  continuous  matrix  of  another  phase,  with 
the  boundary  between  the  phases  given  a  predefined  diffuseness.  The 
second  corresponds  to  a  scattering  from  a  single-phase  two-components 
system.  This  is  to  examine  what  erroneous  a  values  are  likely  to  be 
obtained  when  the  system  actually  contains  no  two  phase  structure 
at  all. 

A.  Structure  containing  spherical  particles 

Construction  of  the  scattering  curves  proceeds  in  three  stages: 

<  J 
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deriving  the  scattering  curve  fox  pin-hole  collimation,  applying 
the  slit  smearing  and  adding  random  statistical  error. 

For  an  ideal  two  phase  system  containing  spherical  particles, 
of  radius  a  distributed  according  to  a  distribution  function  wCa), 
the  scattered  intensity  per  unit  volume  of  the  sample  is 

Iid(s)  =  (Ap)2  f1J<|)2(s,a)w(a)da/J'(4/3)TTa3wCa)da  (17) 

where  f-^  is  the  volume  fraction  of  the  spheres,  and  $  (s,a)  repre¬ 
sents  the  scattering  from  a  sphere  of  radius  a,  given  (Guinier  1963) 
by 

4>(s,a)  =  (4/3)  iTa3U  (2irsa)  (18) 

U(x)  =  (3/x3) (sinx  -  xcosx)  (19) 

For  w(a) ,  we  take  a  triangular  distribution  rising  linearly  from 
aQ-Aa  to  aQ  and  then  falling  linearly  from  a0  to  a0+Aa.  (This  is 
used  here  in  preference  to  the  usual  gaussian  because  it  imposes 
an  upper  and  lower  limit  to  a.)  Eq.  (17)  represents  the  sum  of 
independent  scattering  by  the  spheres  and  ignores  the  effect  of 
inter-particle  interference,  but  this  is  justified  because  in  the 
Porod’s  region  only  the  intra-particle  scattering  makes  recog¬ 
nizable  contribution. 

The  pin-hole  collimation  scattering  curve  I(s)  is  then  ob- 

2 

tained,  according  to  eq.  (1),  by  multiplying  Ij_d(s)  with  H  (s)  to 
account  for  the  diffuse  boundaries  and  by  adding  1^  for  the  'back¬ 
ground.'  Eq.  (4),  valid  for  the  boundary  profile  obtainable  under 

2 

equilibrium  conditions,  is  used  for  H  (s) .  1^  depends  weakly  on  s, 

and  its  value  extrapolated  to  s=0,  for  a  pure  amorphous  material, 
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is  given  CGuinier  1963}  by 

IbCO)  =  p2k.T  C201 

where  p  is  the  electron  density  and  <T  the  isothermal  compress¬ 
ibility. 

The  following  numerical  values  are  used  in  the  computation: 
ip  =  4  x  10-2  electrons/A3 ,  f ^  =  0.25,  aQ  =  1Q0A,  ia  =  20A,  Ib(0)  = 
0.3  electrons^/A^ .  These  are  the  values  which  would  be  found  for 
example  with  styrene-butadiene  block  copolymers.  Three  different 
curves  with  boundary  thickness  t  equal  to  5,  20  and  40A  are  calcu¬ 
lated. 

The  slit  smearing  effect  is  introduced,  by  means  of  ea.  (5) , 
into  the  pin-hole  collimation  intensity  Its)  thus  calculated.  When 
the  infinite-slit  approximation  is  taken,  W^Cu)  is  a  constant,  which 
is  here  set  arbitrarily  equal  to  0.5A.  For  the  'background'  1^,  the 
infinite-slit  approximation  cannot  be  used.  Ib  and  Ib  show  about 
the  same  weak  dependence  on  s,  but  their  relative  magnitudes  depend 
on  the  actual  slit  weighting  function  (u) .  For  the  slit  geometry 
typically  encountered  in  Kratky  cameras,  Ib/Ib  *  0.05  (when  (0)  is 
equated  to  0.5A).  Fmpirically,  the  angular  dependence  of  was 
found  (Vonk  1973,  Todo ,  Hashimoto  and  Kawai  1978,  Hashimoto, 
Shibayama  and  Kawai  1980)  to  obey  the  following  relation  fairly  well 

Ib  =  a  +  bsn  (21) 

with  n  an  even  integer.  The  following  values  of  the  constants  are 
used  in  this  work  for  the  numerical  calculation:  a  =  0.015 

2  ■?  2 

electron  /A“ ,  b  =  600  electron  A  and  n  =  4. 

Finally,  to  simulate  the  counting  statistical  error,  the  exact 
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Its)  values  calculated  at  discrete  s  values  are  replaced  by  random 
numbers  generated  in  accordance  with  gaussian  statistics  of  the 
appropriate  mean  and  variance.  In  the  preset- time  mode  of  data 
collection,  or  in  the  use  of  position-sensitive  detectors,  the 
variance  of  the  error  statistics  is  proportional  to  the  intensity 
at  each  data  point.  The  data  points  in  the  three  constructed  curves 
are  given  statistical  errors  accordingly,  with  their  magnitudes  ad¬ 
justed  to  make  the  standard  deviation  at  the  lowest  intensity  point 
equal  to  1%.  The  resulting  three  curves,  for  t  =  5,  20  and  40A,  are 
designated  SI,  S2  and  S3  for  later  reference,  and  are  shown  in  Fig. 

2. 

B.  Single-phase  structure 

Single  phase  solutions,  containing  polymer  molecules,  can  ex¬ 
hibit  appreciable  electron-density  inhomogeneity  due  to  concentra¬ 
tion  fluctuation  and  give  rise  to  scattering  of  X-ray.  The  intensity 
can  be  especially  pronounced  when  the  structure  is  on  the  verge  of 
phase  transition.  As  a  model  scattering  curve  from  such  a  system, 
we  take  the  theoretical  expression  derived  by  Leibler  (1980)  for  a 
homogeneous  block  copolymer  system  above  the  transition  temperature 
of  domain  formation.  The  intensity  of  scattering  given  by  eg.  (IV-5) 
of  his  work  depends  on  the  parameters  x  (the  interaction  parameter) , 

R  (the  radius  of  gyration) ,  N  (the  number  of  segments  per  molecule) 
and  f  (the  fraction  of  component  1) .  To  express  the  intensity  in 
electron  unit  per  unit  volume,  Leibler 's  expression  has  to  be  multi¬ 
plied  by  (Ap)^Vu  where  Ap  is  the  electron  density  difference  between 
the  components  and  Vu  the  volume  of  a  polymer  segment.  The  follow¬ 
ing  numerical  values  for  these  parameters  are  chosen  for  computation: 
f  =  0.25,  R  =  100A ,  N  =  104,  Vu  =  10  A3  and  Ap  =  0.04  electrons/A3 . 
Two  curves  are  constructed,  one  for  XN  =  17  and  the  other  for  ;<  =  0. 


The  second  corresponds  to  a  very  high  temperature  and  the  first  to 
a  temperature  just  a£ove  the  transition  temperature. 

Slit-smearing  integration  is  then  performed  on  these  two  curves, 
the  'background'  is  added,  and  the  statistical  errors  are  intro¬ 
duced,  all  in  the  same  way  as  described  in  section  A  above.  Fig.  3 

shows  the  two  curves  finally  obtained.  They  will  be  referred  to  as 
0 

LltXN  =  17)  and  L2  (X  =  0). 

XV.  Error  in  the  Evaluated  Boundary  Thickness 

The  method  of  determining  the  boundary  thickness,  described  in 
Section  II,  is  applied  here  to  the  synthesized  scattering  curves  to 
see  how  well  the  original  thickness  values  are  recovered. 

A.  Effect  of  statistical  error 

Fig.  4  shows  the  plot  of  ln[(.I  -  It)  s  *  3  vs.  s  for  the  three 
synthesized  curves  SI,  S2  and  S3.  The  'background'  Ifc  subtracted 
is  the  exact  one  which  was  initially  added  in  the  synthesis,  i.e., 
that  given  by  eq.  (21)  with  a  »  0.015,  b  =  600  and  n  =  4.  The  addi¬ 
tional  error  which  arises  when  a  slightly  erroneous  Iv,  is  subtracted 
is  to  be  examined  later.  Here  the  uncertainty  in  the  value  of  t 
Cor  a)  determined  arises  solely  because  of  the  large  scatter  in  the 
individual  data  points  which  becomes  aggravated  after  the  exact  I^ 
is  subtracted.  The  angular  range  of  data  for  which  the  slope  can 
be  evaluated  is  very  limited,  especially  in  the  case  of.  S3.  The 
data  in  the  low  angle  region  deviates  from  Porod's  law  because  of 
the  interference  effect  reflecting  the  shape  of  the  particles,  'while 
the  high  angle  region  cannot  be  used  because  of  the  excessive  scatter 
of  the  points.  The  arrows  in  Fig.  4  indicate  the  range  of  s  used 
for  the  evaluation  of  the  slope.  The  boundary  thickness  calculated 
from  the  slope  by  ec.  (.14)  is  3 ,  19  and  30  A  for  SI,  S 2  and  S3 
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respectively.  These  values  are  considerably  different  from  the 
expected  values  5,  20,  and  40A.  For  curve  SI,  the  angular  range 
indicated  in  Fig.  1  corresponds  to  X  =  0.7-1. 3  when  s  is  converted 
to  x  by  eq.  (9)  with  a  C=  rrt/4/J)  =3.6.  The  lower  end  of  the  x 
range  thus  extends  slightly  beyond  the  validity  of  the  approximation 
given  by  eq.  (.12)  ,  and  this  partly  explains  why  the  t  value  obtained 
is  too  large.  For  S2  and  S3  the  error  in  t  probably  arises  simply 
from  the  large  uncertainty  in  the  determination  of  the  slope. 

To  see  whether  a  substantially  reduced  statistical  error  in 
the  data  can  lead  to  a  more  accurate  value  of  the  boundary  thick¬ 
ness,  additional  scattering  curves  similar  to  SI,  S2  and  S3  are 
synthesized  with  0.1%,  instead  of  1%,  random  counting  error.  The 
plots  according  to  eq.  (12)  for  these  new  curves  are  given  in  Fig. 

5.  Here  the  ranges  of  s  at  which  the  slope  can  be  evaluated  are 
wider  and  extend  to  larger  angles.  The  evaluated  t  values  are  6, 

19  and  36A,  which  approach  the  true  values  much  more  closely.  Even 
here,  however,  it  would  be  difficult  to  claim  an  accuracy  better 
than  10%  in  the  determination  of  t. 

B.  Effect  of  erroneous  background  subtraction 

In  constructing  the  plots  shown  in  Figs.  4  and  5,  the  same 
'background'  as  originally  added  in  the  synthesis  of  the  scattering 
curves  was  subtracted.  In  practice,  however,  the  'background'  has 
to  be  estimated  from  analysis  of  high  angle  data.  Error  in  the 
estimated  background  can  arise  not  only  because  of  the  scatter  in 
the  data,  but  also  because  the  valid  analytical  form  for  the  back¬ 
ground  intensity  curve  is  usually  not  known.  Several  workers  (Vonk 
1973;  Todo,  Hashiraoto  and  Kawai  1978;  Hashiraoto,  Shibayama  and 
Kawai  1980;  Roe,  Fishkis  and  Chang  1981)  found  the  empirical  eq.  (21) 
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to  fit  the  observed  background  fairly  well,  while  P.athje  and  Ruland 
(.197  61  and  Wiegand  and  Ruland  Cl 97 9)  used  the  exponential  form 

r6  =  a  expCbs2).  (22) 

Among  the  four  possible  analytical  forms,  i.e.,  eq.  (21)  with  n  =  2, 

4  and  6  and  eq.  C221 , 'more  than  one  may  give  apparently  satisfactory 
fit  to  the  data  within  the  error.  The  somewhat  arbitrary  choice 
of  the  analytical  form  can  then  influence  the  value  of  the  boundary 
thickness  finally  evaluated. 

In  the  synthesis  of  curves  SI,  S2  and  S3,  the  background  was 
added  by  eq.  (21)  with  n  =  4.  Fig.  6  gives  the  plot  of  these  scat¬ 
tering  curves  against  s6  as  the  abscissa,  and  shows  that  eq.  (21) 
with  n  =  6  is  a  fairly  satisfactory  approximation.  Here  the  data 
extends  up  to  s  =  0.05.  After  the  backgrounds  evaluated  according 
to  the  straight  lines  drawn  in  Fig.  6  are  subtracted  from  SI,  S2 
and  S3,  the  plots  of  ln[ (I  -  IjO  s  *  ]  against  s  are  constructed  in 
Fig.  7.  The  boundary  thicknesses  evaluated  from  the  slopes  are  12, 

22  and  49A,  respectively  all  substantially  exceeding  the  true  values 
5,  20  and  40A. 

In  Fig.  8  curves  SI,  S2  and  S3  are  replotted  to  show  that  eq. 

(22)  is  also  a  reasonably  good  approximation.  Here  again  s  extends 
up  to  0.05.  After  subtraction  of  the  background  given  by  the  straight 
lines  in  Fig.  8,  the  plots  to  evaluate  the  boundary  thickness  are 
constructed  in  Fig.  9,  from  which  extremely  erroneous  t  values  of  8, 

11  and  9A  are  obtained  for  Si,  S2  and  S3,  respectively. 

There  are  some  latitudes  in  the  values  of  a  and  b  which  are 
determined  from  such  plots  as  given  in  Figs.  6  anc  3.  The  bound¬ 
ary  thickness  values  finally  evaluated  are  less  sensitive  to  the 
uncertainties  in  these  constants,  and  are  found  to  vary  less  than 


10%  in  most  cases.  This  fact  probably  led  to  the  occasional  state¬ 
ments  in  the  literature  that  the  boundary  thickness  values  deter¬ 
mined  were  not  very  sensitive  to  the  method  of  background  subtrac¬ 
tion. 

C.  Illusory  boundary  thickness  in  a  single  phase  structure 

An  experimental  scattering  curve  from  a  specimen  which  does 
not  contain  a  two  phase  structure  may  nevertheless  sometimes  yield 
an  illusory  boundary  thickness  if  the  procedure  is  blindly  applied 
to  it.  The  scattering  curves  synthesized  in  Section  IIIB,  on  the 
basis  of  the  theory  by  Leibler  (1980)  of  a  homogeneous  (single-phase) 
block  copolymer,  are  utilized  here  to  examine  this  possibility. 

Fig.  10  shows  the  two  curves  Ll  and  L2  plotted  against  s4,  to 
evaluate  the  background  in  accordance  with  ec.  (21)  with  n  =  4.  The 
values  of  the  constants  a  and  b  in  eq.  (21)  obtained  from  Fig.  10 
are  0.0196  and  280,  respectively,  which  are  appreciably  different 
from  0.015  and  600  used  in  the  synthesis  of  these  curves  initially. 
The  component  of  scattering  due  to  concentration  fluctuation  decays 
rather  slowly  with  s,  and  some  of  it  is  erroneously  included  in 
the  background  evaluated  in  the  above  manner .  The  erroneous  back¬ 
ground  is  subtracted  from  Ll  and  L2  and  Fig.  11  shows  the  plot  con¬ 
structed  for  the  determination  of  the  boundary  thickness  in  the 
usual  manner.  Quite  possibly  one  could  be  misled  to  taking  the 
slope  as  given  by  the  straight  line  there  and  come  to  the  conclusion 
that  t  is  equal  to  around  10A. 

V.  Discussion 

The  numerical  examples  given  in  the  previous  section  demon¬ 
strate  that  under  favorable  conditions  the  boundary  thickness  t  can 
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be  evaluated  correctly.  They  also  point  out  the  need  to  guard  against 
the  possibility  of  obtaining  seriously  erroneous  results. 

As  discussed  in  the  introduction,  the  basic  problem  is  that  the 
two  types  of  deviations  from  Porod's  law  due  to  the  diffuse  boundary 
and  the  density  fluctuation  produce  opposite  effects,  which  can  be 
separated  from  each  other  only  with  difficulty.  As  the  numerical 
examples  in  Section  IIIB  illustrate,  assumption  of  wrong  analytical 
forms  for  the  background  intensity  can  lead  to  grossly  misleading 
values  of  the  boundary  thickness.  The  danger  of  the  wrong  back¬ 
ground  equation  can  be  alleviated  by  extending  the  collection  of  the 
data  to  higher  angles.  This  may  not  eliminate  it  completely,  how¬ 
ever,  since  there  is  never  a  guarantee  that  an  empirical  equation 
giving  a  good  fit  at  high  angles  is  indeed  a  valid  one  in  the  low 
angle  region  of  importance  for  the  boundary  thickness  determination. 

If  each  of  the  components  forming  the  two  phase  structure  can  be 
obtained  in  pure  form,  more  reliable  background  can  be  obtained  by 
measuring  the  scattering  from  the  pure  materials  separately  and  cal¬ 
culating  the  weighted  average  of  the  two. 

The  statistical  error  of  counting  present  in  the  data  makes 
the  evaluation  of  the  thickness  doubly  difficult.  First,  it  makes 
the  separation  of  the  background  uncertain.  Second,  even  when  the 
background  can  be  subtracted  correctly,  the  scatter  in  the  data 
reduces  the  precision  of  the  obtained  t  value  markedly.  The  ex¬ 
amples  given  in  Section  IIIA  show  that  a  statistical  error  cor¬ 
responding  to  1%  standard  deviation  in  the  background  region  is 
barely  tolerable,  if  an  accuracy  in  t  better  than  20%  is  desired. 

One  should  probably  aim  at  achieving  a  statistical  error  of  less 
than  0.1%  in  the  background  —  an  aim  which  may  not  easily  be 
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attained  in  practice.  To  accumulate  10^  counts  per  channel  in  the 
background  with  amorphous  polymers,  it  may  take  several  days  with 
a  Kratky  camera  fitted  with  a  one-dimensional  position  sensitive 
detector  and  operating  on  a  rotating  anode  tube.  Accumulation  of 
such  high  counts,  for  only  a  few  samples,  may  also  alter  the  detec¬ 
tion  efficiency  characteristics  of  the  position-sensitive  detector 
(JRoe,  Chang,  Fishkis  and  Curro  19^1)  noticably.  An  important  factor 
affecting  the  error  in  t  is  the  relative  magnitudes  of  the  Porod  Law 
intensity  and  the  background  intensity.  Collimation  with  a  long  slit 
will  increase  the  latter  more  than  the  former,  and  therefore  short 
slit  or  pin-hole  collimation  would  be  advantageous  under  otherwise 
comparable  conditions. 

The  conclusions  derived  in  this  work  are  based  on  numerical  re¬ 
sults  obtained  with  synthesized  scattering  curves.  Although  they 
are  constructed  to  model  typical  experimental  curves  in  realistic 
detail,  there  would  of  course  be  experimental  results  which  differ 
greatly  from  them.  In  such  cases  the  assessment  of  the  error  given 
in  this  v*ork  may  not  be  exactly  valid.  One  may  then  be  required  to 
undertake  an  assessment  of  the  error  in  the  evaluated  boundary  thick¬ 
ness  by  repeating  what  is  described  in  this  work,  that  is,  by  syn¬ 
thesizing  a  scattering  curve  embodying  the  structure  under  study  and 
the  slit  geometry  used  and  by  examining  how  well  the  boundary  thick¬ 
ness  built  into  the  scattering  curve  can  be  recovered. 
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Legend  to  Figures 


Fig.  1.  The  plot  illustrates  that  eq.  C12]  approximates  g(x)  well 
for  x  between  0.8  and  6.0.  here  g(x),  given  by  eqs.  CIO) 
and  til) /  is  the  function  representing  the  Porod's  law  in¬ 
tensity  for  structures  having  diffuse  boundary  in  the  case 
of  infinite-slit  geometry. 

Fig.  2.  The  plots  show  the  synthesized  scattering  curves  which  simu¬ 
late  the  SAXS  patterns  obtainable  with  infinite-slit  colli- 
mation  from  samples  containing  spherical  domains  (with  a 
triangular  distribution  of  the  radius  between  80  and  120A) 
imbedded  in  a  continuous  matrix  of  another  phase.  The  back- 

4 

ground  intensity  in  the  form  a+bs  is  assumed.  Curves  SI, 
S2,  and  S3  are  for  the  boundary  thicknesses  of  5,  20  and 
4 0A  respectively.  The  statistical  error  added  corresponds 
to  1%  standard  deviation  at  the  point  of  lowest  intensity. 

Fig.  3.  The  plots  show  the  synthesized  scattering  curves,  which 
simulate  the  SAXS  pattern  obtainable  with  infinite  slit 
collimation  from  a  sample  containing  a  homogeneous,  single¬ 
phase  block  copolymer,  ^he  theoretical  equation  derived  by 

4 

Leibler  (1980)  is  utilized,  but  the  background  (a+bs  ) , 
slit-smearing  and  statistical  errors  are  also  added.  Ll 
is  for  XN  =  17  and  L2  for  XN  =  0. 


Fig.  4.  Curves  SI,  S2  and  S3  are  plotted  in  accordance  with  the 

method  described  in  Section  II  for  evaluation  of  boundary 
thickness.  (For  the  sake  of  legibility,  the  plots  for  S2 
and  S3  are  displaced  downward  by  multiplication  of  0.1  and 
0.01  respectively.)  The  exact  background,  as  was  initially 
added  in  the  construction  of  the  synthesized  curves ,  was 


Fig .  5 . 


Fig.  6. 


F  ig .  7  . 


Fig.  8. 


subtracted  from  them  before  making  the  plot.  The  range  of 
s  at  which  the  slope  is  evaluated  is  indicated  5v  arrows. 

The  thickness  t  obtained  is  8,  19  and  3 CIA  for  Si,  S2  and 
S3,  respectively. 

These  plots  are  similar  to  those  shown  in  Fig.  4,  except 
that  the  curves  Si"*,  S2' ,  and  S3-*  shown  here  are  ones  con¬ 
structed  with  statistical  errors  corresponding  to  0.1%  at 
the  lowest  intensity  data  point.  The  much  smaller  scatter 
enables  the  slope  to  be  evaluated  in  wider  ranges  of  s  ex¬ 
tending  into  larger  angles.  The  thickness  values  evaluated 
6,  19  and  36A  now  agree  better  with  the  true  values. 

The  plots  illustrate  that  a  reasonably  good  approximation  to 
the  background  can  be  obtained  even  when  a  wrong  analytical 
equation  (a+bs6,  instead  of  the  correct  one  a+bs4)  is  assumed. 
The  plots  for  S2  and  S3  are  displaced  downward  for  legibility 
by  subtraction  of  a  constant  term  0.01  and  0.02  respectively. 
The  erroneous  backgrounds,  evaluated  according  to  Fig.  6, 
were  subtracted  from  the  scattering  curves  SI,  S2  and  S3, 
before  constructing  the  plot  to  evaluate  the  boundary  thick¬ 
ness.  (.The  plots  for  S2  and  S3  are  displaced  downward  for 
legibility.)  The  thickness  t,  evaluated  from  the  slopes 
indicated  by  the  straight  lines  drawn,  are  12,  22  and  49A 
for  the  three  curves  respectively,  all  larger  than  the 
true  values. 

The  plots  illustrate  that  the  background,  originally  synthe- 

4 

sized  according  to  a+bs  ,  can  be  fitted  fairly  well  by 
2 

aexp  (is  )  .  (The  plots  for  S2  and  S3  are  displaced  dcwnv'ard 
for  legibility  by  multiplying  by  0.5  and  0.25  respectively.) 
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Pig.  9.  The  erroneous  backgrounds,  evaluated  according  to  Pig.  8, 


were  subtracted  from  the  scattering  curves  SI,  S2  and  S3, 
before  constructing  the  plot  to  evaluate  the  boundary 
thickness.  CThe  plots  for  S2  and  S3  are  displaced  down¬ 
ward  for  legibility 1 .  The  thickness  t,  evaluated  from  the 
slopes  indicated  by  the  straight  lines  drawn,  are  8,  11 
and  9A  respectively,  tbe  latter* two  being  grossly  in  error. 

Fig.  10.  The  single-phase  scattering  curves  Ll  and  L2  were  originally 
synthesized  with  the  background  0.015+600s4.  The  plots 
shown  here  illustrate  that,  despite  the  correct  equation 
with  a  4th  order  term  assumed  for  the  background,  the  best 
fit  obtained  differs  greatly  and  is  given  by  0.020+280s4. 
CThe  plot  for  L2  is  displaced  downward  by  0.01  for  legi¬ 
bility.  ) 

Fig.  11.  The  single-phase  scattering  curves  Ll  and  L2,  from  which 
the  erroneous  backgrounds  evaluated  in  Fig.  10  were  sub¬ 
tracted,  are  plotted  according  to  the  method  to  be  used 
for  evaluation  of  the  boundary  thickness.  CThe  plot  for 
L2  is  displaced  downward  by  multiplication  of  a  factor  0.1.) 
The  downward  trend  of  the  data  points  at  high  angles,  pro¬ 
duced  by  subtraction  of  the  incorrect  background,  can  erro¬ 
neously  be  interpreted  to  mean  the  presence  of  phase  bound¬ 
ary  of  thickness  11A. 
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